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Abstract. A sample of carbon-rich stars (C-stars) in the Small Magellanic Cloud 
(SMC) was selected from the combined 2MASS and DENIS catalogues on the basis 
of their J — Kg colour. This sample was extended to include confirmed C-stars 
from the Rebeirot et al. (1993) spectroscopic atlas. In this combined sample (N = 

1152), a smaller number (N = 1079) were found to have MACHO observations. For 
this sub-sample, light curves were determined and 919 stars were found to have 
high quality light-curves with amplitudes of at least 0.05 mag. Of these stars, 
only 4% have a well-defined single period - most of these have multiple well- 
defined periods, while 15% have highly irregular light-curves. The distribution of 
the logarithm of the period versus magnitude, colour, period ratio (if applicable), 
and amplitude was analyzed and compared with previous works. Variable C-stars 
are distributed in three sequences: B, C, and D from Wood et al. (1999), and do 
not populate sequences with periods shorter than logP ~ 1.5. Stellar ages and 
masses were estimated using stellar evolutionary models. 

Key words. Stars: AGB and post-AGB - Stars: variables: general - (Galaxies:) 
Magellanic Clouds 

1. Introduction 

Asymptotic Giant Branch (AGB) stars can be separated into two classes based on their 
spectra: oxygen-rich (0-rich or M-stars) and carbon-rich stars (C-rich or G-stars). 
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M-stars have more oxygen than carbon in their atmospheres {CjO < 1), while C- 
stars display enri chment [CjO > 1) due to dredge-up caused by thermal pulsation 
l lben fc Renzini 198, IT . These thermal pulses also lead to mass loss, as well as to luminos¬ 
ity variations with periods of '^lOO days or longer and peak-to-peak amplitude variations 
up to a few magnitudes at visual wavelengths. 

All stars are oxygen-rich when they enter the AGB phase. Whether or not they 
become C-stars depends primarily on the efficiency of the third dredge-u p and the extent 


Iben 1981: 


Marigo et al. 1999^ . In metal-poor 


and time-variation of the mass-loss (e.g. 
stars, fewer atoms are necessary to change the envelope from oxygen to carbon 
dominated iCjO > 1); therefore, fewer thermal pulses are needed to convert an M-star 
into a C-star. Conversely, mass loss is expected to be stronger in metal-rich stars, leading 
to shorter ACB and C-star phases. 

In the past, it was thought that both oxygen-rich and carbon-rich Mira 
variables follow a well-determined period-luminosity (PL) relation in the near- 
IR regardless of the host sy stem mean met allicity or type, e.g. in a glob¬ 
ular cl uster ( |FeaslL£t_ah2QQ ^, dwarf galaxy ijClass fc Lloyd Ey ans spira l 

galaxy IjClass et al. IQQfiHyan Leeuwen et al. 199711 . or elliptical galaxy lllieikuba 200# . 
Howeyer, the ayailability of long-term photometric monitoring data proyided 


by the microlensing obserying projec ts, e.g. MAC HO • Alcodc^LaLAMd) , OCLE 


iZebruh et al. 2001: 


Udalski et al. 1997|1 . and EROS IjAubourg et al. 199,11 . and large- 


scale near- infrared (NIR) photo metric suryeys, like the Two Micron All Sky Suryey 
(2MASS, 

(DENIS, 


Skrutskie_et_alM9971) and the Deep Near-Infrared Southern Sky Suryey 
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Enchtein et al. 1997|1 haye opened a new window on this issue and reyealed 


that AGB stars lie on multiple parallel sequences in the PL diagram (Cook et al. 1997; 
Wood et al. 1999). Furthermore, red giant branch (RGB) stars at the tip of the RGB were 
also found to yary (Ita et al. 2002, Kiss & Bedding 2003). These results haye proyided 
new and significant constraints for theoretical pulsation models. 

Differences between 0-rich and C-rich Long-Period Variables (LPVs) haye also been 
found. Using a small sample of AGB stars in the SMC obseryed by the Infrared Space 
Obseryatory (ISO), 2MASS, DENIS and MACHO, Cioni et al. (2003) concluded that the 
period distribution of C-stars peaks at about 280 days. They also noted that C-stars haye 
a larger amplitude with respect to M-stars, contrary to what was de riyed for the LMC 
AGB stars, where both types showed a similar amplitude distribution llCioni et al. 200111 . 
Studying a much larger sample of C and M LPVs in both Magellanic Clouds, Ita et al. 
(2004b) c onfirmed that O- and C-rich Miras follow different period vs. (J — K) colour 
relations IIFeast et al. 1989^ . that C-rich Miras tend to haye greater /-band amplitudes 
at redder J — K colour, and that the amplitudes of 0-rich Miras are independent of 
colour. Groenewegen (2004, G04) reached similar conclusions. 
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However, the studies by Cioni et al. (2001) and Ita et al. (2002) were limited to LPVs 
with P < 1000 days, since the OGLE-II observations only span a time-baseline of about 
1200 days. More recently, Fraser et al. (2005) presented an analysis of the eight year 
light-curve MACHO data for LPVs in the LMC and found that C-stars occupy only two 
of the sequences in the period-luminosity diagram. Furthermore, dust-enshrouded stars 
are located in the high-luminosity ends of the both sequences. 

In this paper, we have extended the work of Cioni et al. (2003) and complemented 
the work of Fraser et al. (2005) by investigating the variability properties of all C-stars 
observed by MACHO in the SMC. We used the Master Catalogue of stars toward the 
Magellanic Clouds (MC^) by Delmotte et al. (2002) and Delmotte (2003) to identify 
C-stars in the SMC. This catalogue provides a cross-correlation between the DENIS 
Catalogue towards the Magellanic Clouds (DCMC - IJKg) and the 2nd Incremental 
Release of the 2MASS point source catalogue {JHKs) covering the same region of the 
sky. C-stars were selected statistically on the basis of their red J — Kg colours {J — Kg > 
1.33; Cioni et al. 2003). This selection was checked through cross-correlation with the 
Rebeirot et al. (1993, hereafter RAW93) catalogue of spectroscopically confirmed C-stars 
ISect. l^ni . C-stars found spectroscopically by RAW93 but with J — Kg < 1.33 were later 
included in the sample. In Sect. [3 we present the selection of C-stars and the extraction 
of the corresponding light-curves from the MACHO database. Section |3 discusses the 
analysis and resulting light-curve parameters. The {logP, Kg), (J — Kg, Kg) and other 
diagrams are discussed in Section while Section [3 concludes this work. Details about 
the method developed to determine periods and amplitudes and the quality assessment 
of the data and of the relevant parameters are given in Appendix A. 

2. Photometric Data and Light-Curves 

2.1. C-stars photometric selection 

The MC^ catalogue, containing the cross-correlation between DENIS and 2MASS sur¬ 
veys, as well as optical UCACl and GSC2.2 catalogues toward the LMC, was published 
by Delmotte et al. (2002). Here, we use its extension to the SMC (Delmotte 2003) and the 
near-IR information of the catalogue only. The region confidently populated by C-stars 
{J — Kg > 1.33 mag and Kg < 12 mag) contains a total of 1657 stars. 805 stars within the 
MACHO fields satisfy the photometric criterion. The MACHO project observed 6 fields 
(each of 0.49 deg^), covering the densely populated bar of the SMC of approximately 3 
square degrees in total. 

Figuredshows the near-IR Kg vs. J — Kg colour-magnitude diagram (CMD) of stars 
from MC^ and within the MAC HO helds. Evolve d AGB stars occupy the region above 


the RGB-tip at Kg ^ 12.6 mag I Cioni et al. 2000 1 and J — Kg ^ 0.97 mag. At Kg ~ 12 
mag the split into two branches is signihcant, though it starts already at Kg ^ 12.5. 
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Fig. 1. Near-IR CMD of SMC stars from 2MASS within the observed MACHO fields 
(dots). The C-stars region is marked by two perpendicular dashed lines: C-stars with 
J — Ks <2 mag are indicated with blue circles; red triangles refer to obscured AGB stars 
(J — Ks > 2 mag). The position of the RGB-tip is also indicated. In the upper right 
corner we show the number of sources with Ks < 12 vs. J — Ks colour; note that the 
branch of C-stars is well separated. See the electronic edition of the Journal for a colour 
version of the figure. 


C-stars populate the well-extended tail toward red colours, while M-stars lie along the 
almost vertical sequence at J — Ks ^ 1.2 mag. They reach a maximum luminosity of 
Ks — 10 mag, except for 3 stars that have Ks — 9 and are likely to be more massive 
0-rich stars (G04). Thus, these 3 have been excluded from the sample. Dust-enshrouded 
AGB stars are located at redder colours {J — Kg ^ 2 mag). These obscured stars can be 
either G-rich or 0-rich, and spectra are needed to distinguish between the two types. We 
include them in our analysis, unless they are explicitly rejected by RAW93 (see discussion 
of contamination by 0-rich stars below). The small box in the upper right corner shows 
the histogram of sources with Kg < 12 vs. J — Kg colour. Indeed, the branch of C- 
stars is well separated from 0-rich stars at about J — Kg = 1.33. Other approximately 
vertical sequences in the main figure at J — Kg < 0.97 are populated by either foreground 
galactic stars or red supergiants and upper main-sequence stars that belong to the SMC 
(i.e. Nikolaev & Weinberg 2000). 

Figured displays the J — H vs. H — Kg colour-colour diagram of the SMC sources 
within the MACHO fields. This diagram is useful for iden tifying stars with large infrared 


excess IjBessell fc Brett lOSS : 


Nikolaev fc Weinberg 200fl|l . The main feature in the dia- 
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Fig. 2. J — H vs. H — Kg diagram of SMC stars within the MACHO fields with Kg < 14 
mag. The dashed line corresponds to J — Kg = 1.33 mag. Symbols are as in Fig. ^ See 
the electronic edition of the Journal for a colour version of the figure. 

gram is an extended branch to the right side of the dashed line that marks our selection 
of C-stars (J — Kg = 1.33 mag). It has been recognized as the locus of TP-AGB stars 
(e.g. Marigo et al. 2003). Circles and triangles, as in Fig. ^ correspond to stars in our 
sample. The less populated region at J — H ^1.2 and H — Kg ^ 0.8 mag corresponds 
to obscured AGB stars. Along this branch a small contamination of stars with Kg > 12 
mag is present (small dots). A handful of them might also be C-stars with very large 
amplitudes caught at minimum light or faint extrinsic C-stars (Westerlund et al. 1995). 



2.2. C-stars confirmed spectroscopically 


Our photometric selection criteria of C-stars were checked against the spectroscopically 
confirmed C-stars in the RAW93 catalogue. Among 1707 stars listed by RAW93 we 
rejected 27 that were classified by the authors as “doubtful” (flag=l in Column 8 of 
their Table 4). Then, we cross-correlated the coordinates of the remaining 1680 C-stars 
with stars in the entire MC^ catalogue brighter than the RGB tip (i.e. Kg = 12.6 mag), by 
adopting a searching radius of 3". The absolute distance between the 2MASS counterpart 
of the RAW93 sources is shown in Fig.|2li- We found that 1275 C-stars in the RAW93 
catalogue have a 2MASS counterpart within 3", and the peak of the distance distribution 
is at c? 1.2". The majority of the spectroscopically confirmed C-stars are well matched 
to s tars in our sample when the same colour and IFs-magnitude criteria are adopted (see 


also 


Cioni et af. 2000|i . Restricting the area to those fields observed by MACHO, there 
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Fig. 3. Panel (a): Absolute distance between 2MASS and RAW93 stars in the MC^ 
catalogue. Panel (b): Kg vs. J — Kg CMD of all SMC stars in the MACHO field (small 
black dots). C-stars included in the RAW93 catalogue and within the MACHO fields 
are plotted as red squares. Blue crosses represent the sample of photometrically selected 
C-stars. Near-IR photometry comes from the 2MASS survey. Panel (c): histogram of 
C-stars in our photometrically selected sample (empty) and in the RAW93 catalogue 
(dashed). See the electronic edition of the Journal for a colour version of the figure. 

are 931 C-stars in common between RAW93 and MC^. The location of these stars in the 
near-IR CMD is presented in Fig. |3). 

We find that 73% of 802 C-stars photometrically selected have confirmed C-type 
spectra. Our photometric selection identifies more C-stars at J — Kg ^ 2 mag with 
respect to RAW93. This might be because 1) some stars with these colours might be O- 
rich AGB stars, and 2) some obscured stars were probably below the detection limit in 
the RAW93 survey. Due to metallicity dependence of the C-star life-time, in a population 
with the metallicity of the SMC we expect only few of our non-RAW93 stars to belong 
to the first category. Many spectroscopically confirmed C-stars have a colour bluer than 
J — Kg = 1.33 mag and a magnitude fainter than Kg = 12 mag, overlapping the region 
where 0-rich stars are also present. These C-stars cannot be disentangled using only a 
photometric selection criterion (see also G04). 
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Figure Et shows the histogram of the number of photometrically selected C-stars 
versus J — Kg colour, together with the number of spectroscopically confirmed C-stars 
by RAW93 versus J — Kg. There is a shift between the two distributions suggesting that 
the spectroscopic identification of C-stars is biased to bluer colours. This is not surprising 
- in a sample that was spectroscopically selected using CN bands near 8000 A, Blanco 
et al. (1980) also found more C-stars with bluer colours than J — Kg = 2. 

We find that within the MACHO fields there are 117 stars with 1.33 < J — Kg < 2.0 
in the MC^ that are not present in the RAW93 catalogue. They amount to 17% of the 
total number of MC^ stars with these colours. Because of the bias discussed above this 
is an upper limit to the number of 0-rich stars contaminating this region of CMD. From 
comparison with the spectroscopic sample of G04, we also expect the contamination to 
be low due to the fact that there are only 3 confirmed 0-rich stars in the SMC with 
(J — Ks)o = 1.30, 1.58, and 2.80, respectively. 

Considering the mis-identifications or missed cross-identifications as a result of an 
automatic association criteria, we estimate a contamination of significantly less than 2.7% 
and 7.9%, respectively, based on Loup et al. (2003). 

2.3. MACHO light-curves 

The light-curves of stars in our sample were extracted from the on-line MACHO cata¬ 
logue^. The MACHO project observed the central body of the SMC simultaneously in 
non-standard (i.e broad) blue (Rm) and red (Rm) bands for roughly 8 years (1992-2000). 
On average there are 800-900 observations per filter for most stars. Since we were inter¬ 
ested in the characteristics of the temporal behavior, we used instrumental magnitudes 
and analyzed the photometric variations in both bandpasses. 

2MASS and MACHO coordinates were cross-correlated using a search radius of 3" and 
the nearest, and reddest star from MACHO catalogue was chosen as the counterpart to 
the 2MASS source. Of the 802 photometrically selected C-stars lying within the MACHO 
fields (Fig.0, 25 stars were detected twice with different identification numbers because 
of overlap between adjacent MACHO fields. However, since we chose the nearest star 
to a 2MASS star as a MACHO counterpart, we avoid double or multiple identifications 
of the same variable. Nevertheless we did check that the quality of the light-curves is 
similar. The histogram in Fig. El shows the distance in arcsec between a 2MASS source 
and the corresponding star in the MACHO database. The distribution peaks at d ^ 076, 
and 65% of the stars are within d <V'. For 51 stars the nearest MACHO counterpart is 
at d > 3". Of these, 20 are within 10"(6 within 4") and 31 within 20". The majority of 
these stars are located at the edges of the MACHO fields and have poorer astrometry 
perhaps due to distortions. Moreover, they have bluer instrumental magnitudes than 


^ http://wwwmacho.mcmaster.ca/ 
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Fig. 4. Locations of C-stars in the observed SMC MACHO fields (6 large squares each 
indicated by its MACHO field number). Small black dots are all photometrically selected 
C-stars in the SMC; blue crosses refer to C-stars with \dMACHO — d 2 MASs\ < 3". 
Small red squares are C-stars from the RAW93 catalogue cross-correlated with the MC^ 
catalogue. See the electronic edition of the Journal for a colour version of the figure. 

other stars in our sample. Therefore they have been excluded from our analysis. Finally, 
the photometrically selected sample contains 751 stars with MACHO light-curves. 

By cross-correlating the RAW93 and MC^ sample with the MACHO data set, we 
found that 328 spectroscopically confirmed C-stars are not included in our photometri¬ 
cally selected sample. Thus, the total number of C-star light-curves analysed in the next 
section is 1079 (751-1-328). 

Tabled is an extract of the full table, available electronically at Centre de Donnes 
astronomiques de Strasbourg (CDS)^, and reports the first 15 lines of the cross-identified 
MACHO and MC^ sources. It contains: MACHO, DCMC, and 2MASS identifier (Cols. 1- 
3); right ascension and declination (in degrees) from the second incremental release of the 
2MASS catalogue (Cols. 4 and 5); the positional difference between 2MASS and MACHO 
coordinates in arcsec (Col. 6) and the RAW93 identification number (if appropriate, 
otherwise 0) (Col. 7). 

3. Analysis of Periods and Amplitudes 

An independent Fourier analysis of the Bm and Rm light-curves was performed to search 
for periodicities in the data. The MACHO time-baseline is about 2700 days, more than 


T-1 r I-1-1-1-1 I-1 1 

207 . 


a 



J_1_ I _1_I_I_I_ I _1_I_I_L 


^ http://cdsweb.u-strasbg.fr/ 
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Fig. 5. Histogram of the absolute distance between 2MASS sources and their MACHO 
counterparts. The dashed line is drawn at 3". 


twice that of OGLE-H database used, for example, by Ita et al. (2004a, 2004b), Kiss 
& Bedding (2004) and G04. Resulting periods above 2600 days (logP = 3.4) should be 
considered less significant. Only photometric measurements that are more accurate than 
0.1 mag are used. The m ethod that extracts the lig ht-curve parameters is based on the 


Lomb-Scargle algorithm IILomb 1976 ; 


Scargle 1982 1 as used by Rejkuba et al. (2003). It 


is described in Appendix A.l. 

Initially 751 light-curves of photometrically selected C-stars were analysed. All the 
light-curve fits and their parameters were inspected visually. Based on this inspection, 
and on parameters returned by the light-curve analysis programmes, three quality flags 
were assigned to each light-curve: flag{l) describes the data quality, flag{2) describes 
the light-curve fit quality, and flag{3) describes the detected periodicities. Table El 
summarizes flag values and their meaning while more details are given in Appendix 
A.2 with examples of light-curves associated to a given flag value. Light-curves with 
flag{l) < 2 produce reliable period determinations, as confirmed by the fact that for 
most of them the Fourier analysis has provided good results {flag{2) < 2). Only for some 
stars (8% and 11%, respectively, in Bm and Rm photometry) classified as good light 
curves {flag{l) < 2), uncertain periodicity {flag{2) = 3) or no periodicity {flag{3) = 0) 
was detected because of highly irregular light-curves. 

The same procedure was applied to the 328 stars common to MC^, RAW93 and 
MACHO but missed by the photometric selection. Aliases were identified from the dia¬ 
gram log P vs. Kg magnitude as those periods that create clear vertical paths (see also 
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Table 1. Cross-identification between the MASTER and MACHO catalogs of our sample 
of C-stars. 


MACHO 

DCMC 

2MASS 

a 

<5 

d 

RAW93 

213.15047.194 

J003516.05-732527.7 

0035158-732527 

8.816091 

-73.424309 

.100 

0 

213.15048.4 

J003521.86-732422.8 

0035216-732422 

8.840220 

-73.406326 

2.400 

0 

213.15054.264 

J003526.42-725835.6 

0035263-725835 

8.859856 

-72.976501 

1.150 

0 

213.15046.550 

J003533.80-733032.6 

0035337-733032 

8.890656 

-73.509033 

.750 

0 

213.15051.6 

J003537.30-730956.4 

0035372-730956 

8.905227 

-73.165588 

.840 

0 

213.15048.8 

J003538.53-732441.3 

0035384-732441 

8.910099 

-73.411438 

1.910 

0 

213.15053.2 

J003547.99-730213.7 

0035479-730213 

8.949644 

-73.037003 

.430 

0 

213.15054.103 

J003552.13-725834.2 

0035520-725834 

8.966831 

-72.976120 

.650 

0 

213.15105.8 

J003603.01-732346.3 

0036029-732345 

9.012422 

-73.396103 

1.210 

0 

213.15106.14 

J003612.60-731711.7 

0036125-731711 

9.052299 

-73.286476 

.360 

0 

213.15105.12 

J003616.80-732133.1 

0036167-732133 

9.069947 

-73.359169 

.660 

0 

213.15108.5 

J003628.84-731144.1 

0036287-731144 

9.119937 

-73.195580 

.170 

1 

213.15108.7 

J003630.13-731033.8 

0036300-731033 

9.125388 

-73.175995 

.380 

0 

213.15109.14 

J003633.23-730548.0 

0036331-730547 

9.138298 

-73.096611 

.320 

3 

213.15106.12 

J003648.03-731830.7 

0036479-731830 

9.199909 

-73.308487 

.320 

6 


Appendix A.l). These correspond to periods equal to 1 and 2 years exactly and were 
removed from our analysis. 

The adopted procedure allowed us to define a semi-automatic algorithm to obtain 
the light-curve parameters and access their quality. It is summarized as follows: the best 
fitting period(s) are obtained from Ea. IA.il the quality of the observations is derived from 
flag{l)', the quality of the period determination is evaluated using the spectral power 
that is closely related with a semi-automatic definition of flag{2) (Appendix A.2). Note 
that this procedure can also be applied to stars of a different type (i.e. M-type stars) in 
the MACHO catalogue or to other measurements of stellar variability in a comparable 
sampling. 

Amplitudes related to the main periodicity of light-curve variations were determined 
both from the sinusoidal fit of each light-curve and from the peak-to-peak magnitude 
difference. A comparison of both determinations is given in Appenxix A.3. 

3.1. Results and statistics 

Table 0 lists the parameters derived from our analysis for a total of 1079 stars. Only the 
first ten lines are shown in this paper, but the complete table is accessible electronically 
via CDS. The table contains: MACHO identifier (Col. 1); the quantities of the period- 
amplitude analysis for Bm (first row) and Rm (second row) light curves: (Cols. 2-9): 
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Table 2. Description of the values of different flags for the complete sample of 1079 
C-stars. 




Data Quality 

BLUE 

RED 

flagil) 

0 

excellent 

306 

803 


1 

good 

478 

70 


2 

fair 

136 

20 


3 

noisy or few data 

115 

150 


4 

no data 

44 

36 



Fit Quality 

BLUE 

RED 

flag{2) 

0 

excellent 

124 

101 


1 

good 

191 

258 


2 

fair 

470 

349 


3 

bad 

127 

168 


4 

underivable 

123 

167 


5 

underivable: flag{l)=4: 

44 

36 


Detected Periodicity 

BLUE 

RED 

flag(3) 

0 

no periodicity 

143 

135 


1 

1 period 

558 

608 


2 

2 periods 

311 

292 


3 

one period 

67 

44 



second period uncertain 




400 

g 200 
+ 

G 
O 

of 0 
0,-200 
-400 
4 

2 
03 

^ 0 
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<1 
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-4 

2 2.5 3 3.5 

log (Pj 

Fig. 6. Comparison between the period and amplitude derived in this work and those 
from Cioni et al. (2003). 
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Amplitudes 


Fig. 7. Comparison between the period and amplitude derived in this work and those 
from Ita et al. (2004b) and G04. 


Table 3. Light curves parameters: X is the mean magnitudes; first (second) row refers 
to i3M(.RM)-hght curve. 


MACHO 

< A > 

Al 


^2 

S2 

Pi 

P2 

Pi 

P2 

■^phot 

flaaW 

/* ag ( 2 ) 

flag (3) 

213 . 15047.194 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.00 

3 

4 

0 


.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.00 

2 

4 

0 

213 . 15048.4 

- 8.139 

.020 

-.077 

.027 

-.020 

1002.284 

199.586 

102.100 

46.200 

.26 

0 

1 

2 


- 9.629 

.060 

.011 

.000 

.000 

1658.224 

.000 

83.200 

.000 

.16 

0 

2 

1 

213 . 15054.264 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.00 

4 

5 

0 


- 5.808 

-.743 

-.917 

-.300 

.776 

1696.516 

359.045 

178.700 

96.400 

1.91 

1 

0 

2 

213 . 15046.550 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.00 

4 

5 

0 


.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.00 

3 
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mean magnitude Xq, Ai, Bi, A 2 , i? 2 , Fi, and P 2 in days, power strength of the first "Pi 
and second period V 2 ; Aphot (Col. 10) flag{l), flag{2), and flag{3) values (Col. 11-13). 
Among the 1079 stars the following have passed the data-quality criterium flag{l) < 2: 
919 Bm and 893 Rm light-curves. 
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In the following discussion, only the good 919 Bm light-curves are used, unless ex¬ 
plicitly stated otherwise. Within this sample 785 stars also have flag{2) < 2 and a 
minimum amplitude of 0.05 mag, thus are all variables. About 4% show a very regular 
variation with only one periodicity. The others appear multiperiodic. Some (7%) show 
a well-defined first period with amplitude variations and a less clear second periodicity, 
while others (59%) clearly show two periods (examples are given in Fig. El- 

The availability of accurate photometry and long-term observations has made the 
distinction between regular and semi-regular (SR) variables more and more difficult 
(Whitelock et al. 1997). The stellar light-curves can be as regular in SR as in Mira class, 
but on average SR variables show smaller amplitudes (Cioni et al. 2003). However, due 
to a difficult and rather subjective classification we only distinguish between two broad 
groups: sources which show a clear single periodicity and sources which are multiperiodic. 
About 15% appear to be irregulars {flag{2) = 3), with no clear period. 

Table 0 has 131 stars are in common with Cioni et al. (2003). Fig. Et shows the 
comparison between the periods derived in the present paper and those in Cioni et 
al. (2003). Amplitudes are compared in Fig. Et. We plotted stars with flag{2) < 2 in 
our analysis and stars with Flag < 9 and Flag 5 in the table by Cioni et al. (2003). 
There are 78 stars in common. The mean period and amplitude differences are: Ppresent — 
Pcioni 03 = -14 ± 76 and Apresent “ Acionios = “0.9 ± 0.8. The periods agree within the 
uncertainty in the period determination that is on the order of 5%, except for few stars 
for which Cioni et al. (2003) derive only one period and flag the LPV as multiperiodic, 
while we find 2 periods that are typically shorter. It is possible that additional longer 
periods are present as well. In contrast, amplitudes are systematically different. Note 
that Cioni et al. (2003) define amplitude as the difference between the minimum and 
maximum value of MACHO photometry, which is different from our definition of Aphot 
(see Appendix A.3). This could explain why Cioni et al. amplitudes are systematically 
larger than ours. 

Figure [3 shows the comparison between periods and amplitudes as derived in the 
present paper and those by Ita et al. (2004b) and G04. In both panels, for all the three 
works only stars satisfying our photometric criteria are reported. In Fig. Eh in the case 
of multi-periodic variables the period of G04 corresponding to the largest amplitude and 
our first period are reported. Ita et al. (2004b) only give the predominant period and do 
not analyze multi-periodic light-curves. The present results are in good agreement with 
those by G04, predicting periods as long as 2400 d, while Ita et al. (2004b) do not find 
periods longer than about 690 d. 

In Fig. Eh the OGLE /-band amplitudes of Ita et al. (2004b) and G04 are compared 
with the present i?M“band amplitudes as derived from the light-curves { Aphot ( Rm ), 
see Appendix A.l). In the case of G04 we plot the largest amplitudes from his Table 
2. Although a direct comparison between Rm- and /-band amplitudes is difficult, the 
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log (P) 

Fig. 8. { Ks)q vs . logP relation for all C-stars in the sample (small black dots). Large 
dots are first periods (Pi) for stars with flag{2) < 2 and flag{3) = 1; blue open squares 
Pi for stars with flag{2) < 2 and flag{3) = 2 and 3; red crosses correspond to the 
second period (P 2 ) of stars indicated by squares. The {Ks)o magnitude is the average of 
DENIS and 2MASS measurements, dereddened as explained in the text. Letters identify 
the parallel sequences as defined by Wood et al (1999). 

shape of the histogram illustrating our results is similar to that by Ita et al. (2004b), 
even if we have a larger number of stars with amplitudes smaller than ^0.4. G04 found 
that the majority of C-stars in his sample have A/ < 0.2, while Ita et al. (2004b) have 
almost no stars in the same amplitude range. 

Both Fig. 1^ and Fig. m indicate that the various methods for estimating pulsation 
amplitudes can give systematically different results (see also Appedix A.3). Ita et al. 
(2004b) and Cioni et al. (2003) derive pulsation amplitudes as AX = Xjnax—Xmin, where 
X is the photometric band used, while G04 estimate amplitudes from the sinusoidal fit. 

4. Discussion 

4.1. Log(P) vs. K diagram 

In the past five years new results for the log(P) vs. AT-band mag diagram have been 
obtained as a by-product of microlensing projects such as MACHO, OGLE, and EROS. 
Wood et al. (1999) found that the bright red giant variables in the LMG form four 
sequences: three are the result of different pulsation modes and one, at the longest periods 
(seq. D or the long secondary period sequence - LSP), remains unexplained. Ita et al. 
(2004b) showed that some sequences possibly split into sub-sequences at the discontinuity 
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Fig. 9. Histogram of all the periods found for Bm light-curves, though for Rm light- 
curves results are very similar. The empty histogram shows the distribution of Pi and 
the shaded histogram of P 2 . Letters indicate sequences as in Wood et al. (1999). 

around the tip of the RGB. G04 analyzed spectroscopically confirmed M- and G-type 
stars and concluded that the LSP sequence is independent of evolutionary and chemical 
(G-rich or 0-rich) effects. A comparison of the SMG and the LMC sequences can also 
be found in Gioni (2003). More recently, Schultheis, Glass, and Cioni (2004) compared 
variable stars in the different sequences between the Magellanic Glouds and the NGG6522 
field in the Galaxy. These authors conclude that all three fields contain similar types of 
variables, but the proportion of stars that vary decreases at lower metallicities and the 
minimum period associated with a given amplitude gets longer. 

FigurelHlshows the distribution of the mean {Ks)o mag vs. the logarithm of the period 
obtained by analyzing Bm light-curves. However, the discussion that follows applies also 
to a similar analysis of Rm light-curves. In fact as expected from the comparison shown 
in Fig. IA.6I there are no differences in the periods obtained from the two channels. 

In Fig. IHl all sources in the sample are indicated by small black dots. In addition, 
stars for which only one reliable period was detected {flag{2) < 2 and flag{3) = 1) are 
plotted as larger dots, while for stars with two reliable periodicities detected {flag{2) < 2 
and flag(3) = 2 and 3), the first and the second periods are plotted with blue open 
squares and red crosses, respectively. The ATg-band photometry is the mean of DGMG 
and 2MASS measurements ((ATg)). DGMG magnitudes are corrected for the shift in 
the absolute calibration according to Delmotte et al. (2002). The ATg-magnitude is also 
dereddened according to a SMG mean reddening of E{B — V) = 0.065 ± 0.05 obtained 
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by averaging different measurements llWesterlund 1997^ . By adopting the extinction law 


of Glass (1999) and Rv = 3.1 we obtain Ak^ = 0.02 and Aj = 0.05. 

All periods range in the interval 1.5 ^ log(P) ^ 3.5 (Fig. |^. There are 3 well- 
defined parallel sequences in the PL diagram with a small number of stars lying between 
the sequences. Part of the scatter is due to the fact that the ATs-band is a mean of 
only 2 measurements and part is due, probably, to depth effects in the SMC. The line 
of sight depth of the SMC is estimated to range between 5-20 kpc by various authors 
I Westerlund 199711 . Hence, assuming the SMC distance modulus of (m-M)=19 mag and 
the full depth of 10 kpc, we expect a scatter on the order of ^ 0.4 mag around the 
average. 

Sequences can be identified with B, C, and D from Wood et al. (1999). It is interesting 
to note that C-stars do not populate shorter sequences (i.e. sequence A), a result already 
visible from Fig. 4 of Ita et al. (2004b), and noted in the LMC by Fraser et al. (2005). 
Since the baseline explored here is longer than that of previous works (except for Fraser et 
al. 2005 who analysed MACHO data in the LMC), we find a well populated D sequence, 
about 34% of the variable stars in our sample have a first or second period longer than 630 
days. This should be com pared with previou s results which found (i) 25% of all variable 
AGB stars in the LMC llWood et al. 199911 . but with a very small fraction belonging 
to C-rich LPVs (Fraser et al. 2005), and {ii) 24.6% of all the spectroscopically selected 
C-stars with periods from OGLE-H photometry in the SMC on the sequence D (G04). 
The bias against detection of variables with periods in excess of ~ 800 days in the latter 
work may be why we find more C-stars on this sequence. This shows that very long-term 
monitoring is essential. Even in the case of our sample from MACHO with an 8-year 
time-baseline, it is not clear if the drop in the period distribution at log P ~ 3.3 is real 
or an artifact due to incompleteness at longest periods (see Fig.^J. 


Sequence D is much broader tha n the others, and the nature of its stars is still a 
matter of debate llWood et al. 200411 . Only a few of them, those with J — Kg > 2 (see 
Sect. 14.511 are probably dust-enshrouded AGB stars that could have either carbon or 
oxygen-dominated chemistry, but their number is very small. Thus, as already concluded 
by Wood et al. (2004) from the similarity of the colour variations associated with both 
the primary and secondary periods, dust is unlikely to cause the LSPs. Clearly, a long¬ 
term spectroscopic and photometric monitoring of these stars is necessary to gain some 
insight into the nature of their variability. 


In Fig. Owe compare the histogram of the first detected periods, which represent the 
dominant periodicity in a given star, with the second best period. First periods mainly 
occupy sequences C and D with the peaks at log P = 2.45 and log P = 2.95, respectively. 
Second periods mainly populate sequence B and peak at log P = 2 — 2.3. A weak peak 
is also seen at log P = 3.2, which is clearly part of sequence D. 
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4.2. Log(P) vs. period ratio 

According to Lattanzio & Wood (2003) a star that evolves up the AGB first pulsates 
at low amplitude on sequence A, then with further evolution the sequence B mode will 
become unstable and its amplitude increases, while the sequence A mode amplitude de¬ 
creases. The star will be a multi-mode pulsator having periods in sequences A and B. The 
pulsation amplitude of each mode increases with increasing stellar radius. Subsequently 
a fundamental mode in sequence C will also become unstable. At this moment up to 
three different periods of pulsation can be detected in a given star. However, the com¬ 
peting growth rate of the amplitude of each pulsating mode may shade the detectability 
of a given periodicity. For example for a ^Mq star the fundamental period dominates at 
Mboi > —4.5, while at Mboi < —3.5 the first and second overtone pulsation modes are 
stronger. The star will finally end up as a dust-enshrouded, large amplitude fundamental 
mode pulsator prior to ejection of all its envelope and the beginning of the Planetary 
Nebulae phase. 

For stars with detected multiperiodicity we plot the ratio between the longer to the 
shorter period as a function of the longer one in Fig. cni A well-defined group of stars 
have ratios ranging from 1 to 2. These are shown in zoom-in in the upper left corner of 
the figure. They belong to the sequence C (longer) and B (shorter) periods. Their ratio 
distribution agrees with the scenario proposed by Lattanzio & Wood (2003) for a star 
pulsating in the fundamental, first and second overtone (see their Fig. 53). An extended 
tail up to a ratio of 20 is populated by stars with the longer period on sequence D and the 
shorter period along either sequence C or B. There are no theoretical models available 
at present to explain these high ratios which involve an intrinsic stellar pulsation; thus a 
different nature for long-term modulation needs to be invoked. It is interesting to note 
that most of the sources with Pi < P 2 have lower period ratios and slightly larger values 
of longer periods with respect to sources with Pi > P 2 . 

4.3. Log(P) vs. amplitude 

Fjgiire rm shows the period distribution (of only the first periods Pi ) as a function of am¬ 
plitude in the MACHO bands. Most variables have amplitudes in the optical bands below 
about 0.5 mag. Stars that occupy sequences C and D are clearly separated. Sequence B 
is not present in these diagrams because it is mostly populated by secondary periods. 
It should be also noted that amplitudes belonging to secondary periods are typically 
smaller, and thus this figure should be compared with those of other authors (e.g. G04, 
Ita et al. 2004b) with caution. In addition, the definition of amplitude is not always 
a trivial issue with these highly variable stars that often have variable amplitudes as 
well. We discuss this in Appendix A.3. Here we use amplitudes determined directly from 
photometry (Aphot) as defined in Appendix A.3. 
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Fig. 10. The ratio between the longer to the shorter period versus the longer period. 
Black-filled circles represent stars with Pi > P 2 , while red open squares sources with 
Pi < P 2 - In the upper left corner we show an enlargement of the lowest ratios. 


Both C and D sequences present distribution tails to large amplitude (from a few 
tenth up to 3 mag) values that probably correspond to variables of Mira type. While 
Miras typically occupy sequence C, a few of them can be found on sequence D because 
they have fainter magnitudes due to dust obscuration. 

Figure IT^ shows the period distribution for small (0.05 < AphotiBut) < 0.2; open his¬ 
togram) and large {AphotiBu) > 0.2; shaded histogram) amplitude variables. The two 
distributions are different. Small amplitude variables have two peaks that correspond to 
the period-magnitude relations C and D. Large amplitude variables are more homoge¬ 
neously distributed between log{Pi) = 1.6 and 3.1. Fewer of them have longer periods, 
though at log P ^ 3.1 the sample might be incomplete. 

In Schultheis, Glass & Cioni (2004) the period of both small and large amplitude 
variables defined as in Fig. 1121 increases progressively with decreasing metallicity, even 
though the general period distribution of the two classes is fairly similar. In the present 
paper stars located in sequence C (1.9 < log{Pi) < 2.7) span the full range of periods 
if they are either small or large amplitude variables. Therefore there seems to be no 
indication of a difference in metallicity between the small and large amplitude LPVs 
on this sequence. On the other hand, the majority of large amplitude variables that 
occupy sequence D (2.7 < log{Pi) < 3.5) have log{Pi) < 3.1. A longer monitoring time- 
baseline is necessary to discern whether this lack of longer period large amplitude LPVs 
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Fig.11. Period distribution as a function of amplitude in the two MACHO photometric 
bands. 



Fig. 12. Period distribution {log{Pi)) of stars with 0.05 < Aphot{BM) < 0.2 (empty) 
and Aphot{BM) > 0.2 (shaded). 


on sequence D is real or due to incompleteness, and thus if it could represent shortage 
of lower metallicity stars among the large amplitude variables. 
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Fig. 13. (J — Ks)o vs. the i?M-band amplitude (a) and vs. log (Pi) (6) for all C-stars. 



4.4. Log(P) vs. J — Kg colour 

Figure shows that redder variables have a larger amplitude. In fact Ita et al. (2004b) 
also found that the C-rich regular pulsators (Miras) have larger amplitude the redder 
the star, while the 0-rich Miras seem to have arbitrary amplitudes as a function of 
J — Kg colour. When the J — Kg colour gets redder, periods increase for variables in the 
C sequence (see Figure Et)- This is apparently not the case for stars populating the D 
sequence, which again indicates a different mechanism is responsible for the light-curve 
variations. 

4.5. Properties of stars on the BCD sequences 

We found no difference in the location of stars belonging to the BCD sequences in the 
Kg vs. J — Kg diagram. In Fig. 11 4l we investigated the Kg luminosity function and J — Kg 
colour distribution of stars in each sequence for both the first (empty histogram) and 
second (shaded histogram) detected periods. Most of the stars with J — Kg >2 belong 
to sequence C, and they have only one detected and dominant periodicity. A handful 
of them show a long second period in sequence D. The faintest stars analyzed in this 
sample (12 < Kg < 12.6) have their first period either in sequence B or D and eventually 
a second period in sequence C, which is often the case for stars in sequence D that are 
multi-periodic, while stars in sequence B usually have only one period detected. The bulk 
of the C-star population has a first period that occupies sequence C or D and somewhat 
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Fig. 14. Distribution of the number of sources populating each BCD sequence versus 
it's magnitude and versus J — Kg colour. The empty histogram indicates the first period 
(Pi) and the shaded histogram the second period (P 2 ). 



less sequence P, which is instead largely populated by the second period of these same 
stars. 


4.6. Log(P) vs. Mboi 


FigureElshows the PL-relation for all stars in the sample. Overplotted are the theoretical 
models by Vassiliadis & Wood (1993) for different masses and a mean SMC metallicity 
oi Z = 0.004. To transform Kg magnitudes into M^oi we used the relations by Bergeat, 
Knapik, & Rutily (2002). Their Fig. 1 shows that the bolometric correction is well-defined 
for C-stars with J — Kg ^ 2.1, while for increasing colours the uncertainty becomes larger. 
In our sample of 1079 stars, only 30 have (J — Kg)^ > 2.1 and flag{2) < 2, these sources 
are emphasized in Fig. 1151 as green triangles. These stars have large amplitudes (see 
Fig.HSll which is likely to influence their location in Fig. d 


A comparison of the PL diagram with the models in Vassiliadis & Wood indicates 
that the bulk of C-rich LPVs in the SMC have lower masses than 2.5-3 M^t^, which at 


the metallicity of the SMC corresponds to ages of 0.5-0.3 Gyr llPietrinferni et al. 200411 . 


This agrees with the SFH derived by Harris & Zaritsky (2004), who found an active star 
formation during the past 3 Gyr with three enhanced episodes at 2.5, 0.4, and 0.06 Gyr. 
Few carbon stars show higher masses up to 5 M© corresponding to an age of ^^0.1 Gyr. 
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Fig. 15. Mjjoi vs. logP relation for all stars in the sample and overplotted theoretical 
tracks by Vassiliadis & Wood (1993). Green triangles indicate stars with {J — Ks)o > 2.1 
and flag{2) < 2, while other symbols are as in Fig. |H1 Masses are indicated on the right 
hand side of each track. See the electronic edition of the Journal for a colour version of 
the figure. 

4.7. Spatial distribution 

The spatial distribution of stars in each sequence is investigated. There seems to be no 
indication of a spatial correlation between stars that occupy one or the other sequence 
in the period-magnitude diagram with their location in the galaxy. We also checked if 
there is any difference in the spatial distribution of stars where different masses were 
expected. Again, we do not find clear indication of a different distribution of sources 
with a different initial mass. There are overall less sources with high mass (or young), 
and a lack of these sources in the Northeast, Northwest and Southwest regions compared 
to other locations in the galaxy. The distribution of sources with intermediate and low 
mass is fairly similar. 

5. Summary and conclusions 

This work analyzes and discusses the MACHO light-curves of 1079 C-stars. The sam¬ 
ple consists of 751 photometrically selected stars from the MC^ catalogue according to 
J — Kg > 1.33 mag and AT < 12 mag criteria, and 328 spectroscopically confirmed C-stars 
from RAW93 catalogue. Many of the photometrically selected sample also have C-type 
spectra (RAW93), meaning that for 18% of all the sample we do not know their spectral 
type. However, given the low metallicty of the SMC and a very low number of spectro- 
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scopically known 0-rich AGB stars in the SMC (G04), we expect the contamination by 
0-rich AGB stars to be negligible. For all the stars, we performed Fourier analysis of 
their light-curves identifying up to 2 significant periodicities. 

All the stars for which good quality light-curves exist were found to vary with am¬ 
plitudes of at least 0.05 mag in the two MACHO photometric bands. The analysis was 
carried independently for Bm and Rm light-curves, and the derived periods are identical 
(to within the errors). After a selection based on the quality of the light-curve, signif¬ 
icance of the derived periods, and quality of the periodicity fits, a total of 919 C-stars 
were used in further analysis. 


Carbon stars occupy bright parts of sequences B, C, and D in the Kg —log P diagram. 
None of the stars in our sample have shorter periods characteristic of sequence A, which 
is in agreement with recent studies by Fraser et al. (2005) in the LMC, and Ita et al. 
(2004a,b), and G04 who observed this in the SMC as well. The large majority of the stars 
have their primary (dominant) periodicity on sequence C or D, while more than 2/3 of 
sequence B is populated by secondary periods of those stars whose primary period is on 
sequences C or D. The stars whose primary period is on sequence B are preferentially 
fainter and bluer. This is in agreement with the models that predict change of pulsation 
mode for the LPVs f rom higher overtones tow ards a fundamental mode as they evolve 
along the AGB (e.g. 


Lattanzio fc Wood 2003^ . 


Most of the stars with J — Kg > 2 belong to sequence C, and only a handful of 
these reddest variables show a long second period in sequence D. This is a clue that dust 
obscuration cannot be a cause of the long secondary periods (see IWood et al. 2004 for 


a detailed discussion). The luminosity functions of the three sequences span a similar 
range of magnitudes, but the faint distribution tail of sequence D is more populated. 
Again, we cannot explain it with the dust obscuration, as there are very few stars redder 
than J — Kg > 2 on this sequence. The width of this sequence in the period-magnitude 
diagram is larger than that of other sequences, especially at brighter magnitudes, but 
the lack of stars with a period longer than (log P ~ 3.5) may be due to incomplete time 
coverage. 


Stars belonging to different pulsational sequences are homogenously mixed over the 
SMC area. Their masses were derived from a comparison with the theoretical tracks of 
Vassiliadis & Wood (1993), and the majority of them are indicative of the major star 
formation episode that took place ~ 0.3 — 0.5 Gyr ago. This is in excellent agreement 
with Harris & Zaritsky (2004), who found enhanced star formation episodes at 2.5, 0.4, 
and 0.06 Gyr in the SMC. However, our sample indicates a considerably weaker star 
formation event at younger ages (~ 0.1 Gyr ago), which is expected as these younger, 
and thus more massive, stars are preferentially oxygen-rich, hence not in our sample. 
There is no clear difference in the spatial distribution of the stars with different masses. 
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The very long-time baseline of MACHO observations has allowed us to confidently 
extract long periods up to {log P ~ 3.5). We found that about 10% of the variables 
fall on sequence B, 30% on C, and 34% on D. The latter percentage is higher than the 
25% derived by Wood et al. (1999) and than the 21% derived by G04 in the LMC. It is 
also important to note that in the LMC, Fraser et al. (2005) hnd only a small fraction 
of probable C-stars with periods along sequence D. Monitoring of these variables over 
more than an 8-year period photometrically and spectroscopically is desirable in order to 
discern their nature. According to Wood et al. (2004) these stars belong to the only class 
of bright large amplitude variables whose properties cannot be explained with theoretical 
models at present. 

Appendix A: Determination of Periods and Amplitudes 

A.l. Method 

The method used to determined the first period is already discussed in Rejkuba et al. 
(2003). To detect the second period we used the following fitting function: 



(A.l) 


+ A, 


0 


where the first guess for the second period comes from the second, or sometimes the third, 
strongest peak in the power spectrum. Usually, the hrst period is obvious, while secondary 
peaks are less clear. We noted that semiregular and irregular light curves show several 
peaks in their Fourier spectrum with similar strength. Usually we followed the strength 
sequences of the power values to define the second period. However, the amplitudes of the 
peak are just a starting point. In fact we always inspected the relative results by eye, in 
order to check the reliability of the second assigned period. In some cases we considered 
the third peak more appropriate than the second, which could be instead an alias (a false 
period which seemingly fits the data, as well as the correct period). The alias frequencies 
Va are related to the true frequency vt by Va = ^ where At is the spacing in 

days between the measurements. The pattern of the aliases will therefore depend on the 
true period and the particular At for the star in question. 

We call the first period (Pi) the best fitting period obtained from the first term of 
Eq.lO and Pi the associated value of the power. We call the second period (P 2 ) the 
best fitting period obtained from Eq. EH and P 2 the associated value of the power. 
Amplitudes (Afu) associated with the first period are defined from the sinusoidal ht as: 



(A.2) 
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Fig. A.l. Examples of light-curves with different values of the data-quality flag{l). 
Data with photometric errors 0.05 < ax <0.1 are plotted as black dots; data with 
0.01 < ax < 0.05 are plotted as red squares; and data with ax < 0.01 as green crosses. 
Left (right) panels illustrate Bm {Rm) photometry. See the electronic edition of the 
Journal for a colour version of the figure. 


A.2. Quality assessment 

Table 121 summarizes the values of the three quality flags associated to each light-curve 
and the number of light-curves with a given value. In particular, flag{l) is related to 
the accuracy of the photometric measurements. We defined no.os and no.oi as the ratio 
between the number of observations with, respectively, magnitude error ax < 0.05 and 
< 0.01 and the total number of measurements in each light-curve (i.e. with magnitude 
errors A < 0.1). Thus, values to flag{l) are assigned as follows: 


flag{\) = 0 if no.oi > 0.7 {excellent) 
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Fig. A. 2. Examples of light-curves with different values of flag{2). Left (right) panels 
illustrate Bm (Rm) photometry. For each band we show the periodogram and the light- 
curve in the phase and time domain. See the electronic edition of the Journal for a colour 
version of the figure. 
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{fair)- 

flag{\) = 3 if no .05 < 0.5 and ng 01 < 0.05 
{few or noisy data) 

A sample of light-curves of different types, together with the associated flag{l), is 
shown in Fig. ro The algorithm works well; only for a handful of cases the inspection 
by eye has revealed that data are sparse over several magnitudes even if they have very 
small magnitude error. In these cases we re-assigned flag{l) = 3. The majority of light- 
curves classihed as flag{l) = 2 look like 213.15624.20 i?M-light curve in Fig. IA.l| Only 
in a few cases are very regular and large amplitude LPVs classihed with flag{l) = 2 or 
even 3 (see for example 213.15332.2889, Bm, and Rm light-curve in the Fig. roi . 

The ht quality flag {flag{2)) and the periodicity flag {flag{3)-, see TableEJ are less 
objective and much more related to an inspection by eye than flag{l). In order to assign 
values to these hags, hrst we carefully eye-inspected the hrst period in the phase and 
time diagrams. Then, we looked at the combined periodicities (Eq. Ell). On the basis 
of these inspections we assigned values for flag{2) and flag{3) according to Table|21 In 
Fig. IA.2l we plot light-curves with different values of flag{2) referring only to their hrst 
periodicity. Fig. El illustrates examples of light-curves with diherent values of flag{3). 
In both hgures for each star we plot the corresponding power spectrum and the light 
curve in phase and time domain and overplot the best htting sinusoidal function. 



0 100 200 300 400 500 

Power (P,) 


Fig. A. 4. Histogram of the power of the hrst period (i?M-light curves) as a function of 
the dehned ht-quality hag. The relative number of stars is indicated in each panel. 
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1.5 2 2.5 3 3.5 1.5 2 2.5 3 3.5 


log (Pb) 

Fig. A.5. Distributions of the first period (SM-hght-curves) as a function of the power 
value (left panels) and flag{2) (right panels). 

In order to provide an automatic classification of the fit quality, we analyzed the 
correlation between flag{2) and the strength of the power of the first period (Vi) as 
follows. In Fig. IA.4I the histogram of Vi is reported for different values of flag{2). It is 
clear that starting from flag(2) = 0 to flag{2) = 4 the peak of the distribution moves 
toward lower power values. In particular, for flag{2) = 0 and 1 (Fig. lA.dh l the only 
region populated is that with Vi > 100; the bulk of variables with flag{2) = 2 have 
Vi ^ 100 with only a small number of stars with Vi < 100. The distribution of variables 
with flag{2) = 3 peaks in the range 80 < < 100, and those for which no satisfactory 

period could be derived mainly have Vi < 40. 

Figure ro emphasizes the good correlation between these two quantities. In panels 
(a, h) we plot log{P) vs. Kg (see next section for a discussion) for stars with 60 < < 100 

(a) and stars with Pi > 100 (b). Panels (c, d) show the same quantities but for stars with 
2flag{2) = 3 (c), and flag{2) < 2 (d). We conclude that the selection described in 
panels (6, d) outlines three sequences and the distribution of stars in both panels are 
very similar. These results are related to Bm light-curves, however an inspection of Rm 
light-curves gives similar results. Fig. IA.6I shows an almost perfect correlation between 
the first period obtained from Bm and Rm light-curves. Only stars with flag{2) < 2 
are plotted. For 94% of the objects the residual dispersion of the differences is lower 
than 4%. Only a few points scatter away from the 1:1 line, implying that the primary 
period is the same in both bands. Most of these points occur because the first period 
(here Vi) found from the Bm light-curve corresponds to the second period {V 2 ) found 
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Fig. A.6. Relation between the first periods derived from Bm and Rm light-curves. 

from the Rm light-curve, and confronting only Vi between both light-curves produces 
a point that does not follow the correlation. A few other points occur when the period 
Vi derived from the Bm light-curve is twice that derived from the Rm light-curve. 

Summarizing, an automatic selection criterion (i.e. by using Vi) can be applied with 
a good level of confidence for the analysis of the primary periods of a large sample of 
variables. 

To define flag{3) in a similar automatic way as flag{2), we considered the values 
of the reduced of the sinusoidal fit. In most cases the value decreases when two 
periodicities are considered, even in cases we classified with flag{3) = 1, i.e. when only 
one period could be reliably derived. Perhaps this indicates the complexity of stellar 
pulsation. For these stars flag{3) = 1 does not mean that the light-curve is clearly 
characterized by one periodicity only; rather it implies that the second period is much 
less regular and could not be fitted well with a set of sinusoid functions. 

A.3. Comparison between photometric amplitude and amplitude of the sinusoidal fit 

Figure EH shows the relation between the pulsation amplitude derived from the sinu¬ 
soidal fit {Afit) and that from the peak-to-peak magnitude difference (Aphot) for Bm and 
Rm light-curves separately. The latter value is defined as Aphot = BMirnax) —BM(niiri), 
where BMimax) and BMirnin) are the maximum and the minimum values averaged over 
a few data points in order to avoid spurious detections. For sources of regular period¬ 
icity with excellent observational data (i.e. low photometric errors and flag{l) = 0) 
that clearly show only one periodicity, we discarded the first 5 measurements and took 
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Fig. A.7. Panel (a, b): The residual {Aphot— Afit) is plotted against Aphot- Panel (c, d): 
Histogram of the amplitude is reported for Bm (c) and Rm {d) magnitudes. Panel (e): 
Difference between amplitude from Rm- and i?M-light curves. Only stars with flag{2) < 
3 are shown. 

the mean over the brighter (fainter) ten measurements (see for example the light-curve 
212.15910.2115 in Fig.fOll. Since the majority of sources show less regular light-curves 
with bumps and multiperiodicity, the first brighter (fainter) 150 measurements, which 
could be disturbed by the secondary periodicity with P ~ 150 days, are excluded. Then, 
we computed the “maximum” (“minimum”) value as the average over the next brighter 
(fainter) 40 photometric measurements. 

The sinusoidal fit clearly predicts smaller amplitudes compared to Aphot- 
Quantitatively, AphotiBn) - Afu^BM) ^ 0.12 ± 0.22 and AphotiRm) - AfitiRn) ^ 
0.08 ± 0.25. This is expected in case of deviations from pure sinusoidal variability, due to 
the presence of irregularities in the intrinsic light variation and photometric uncertainties 
caused by measurement errors and variations in seeing. The latter may lead to blending 
and spurious measurements due to the presence of cosmic rays and defects on the CCD. 
As can be seen from Fig. ro this applies to all but the most regular variables with 
flag{2) = 0. 

The histograms of amplitudes in the two MACHO photometric bands are reported in 
Fig. lA.7h . d. The bulk of sources have 0.2 < AphotiBu) < 0.5 and 0.1 < Aphot{BM) < 0.4. 
As expected, amplitudes in Bm are slightly higher than those in Rm (see Fig. lA.Tb b 
The difference in the mean amplitude of the two bands is Aphot{BM) — AphotiRm) = 
0.15 ±0.18. 
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